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We present a measurement of the time-dependent CP-violating asymmetries in decays of neutral 
B mesons to the final states D*^^, using approximately 232 million BB events recorded by the 
BABAR experiment at the PEP-II e + e~ storage ring. Events containing these decays are selected 
with a partial reconstruction technique, in which only the high-momentum n ± from the B decay and 
the low-momentum 7r T from the D* T decay are used. We measure the parameters related to 2/3 + 7 
to be a D * v = -0.034 ± 0.014 ± 0.009 and c l D ,^ = -0.019 ± 0.022 ± 0.013. With some theoretical 
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assumptions, we interpret our results in terms of the lower limits | sin(2/3 + 7)| > 0.62 (0.35) at 68% 
(90%) confidence level. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) quark- 
mixing matrix [1] provides an explanation of CP violation 
and is under experimental investigation aimed at con- 
straining its parameters. A crucial part of this program is 
the measurement of the angle 7 = arg (-V ud V* b /V cd V* b ) 
of the unitarity triangle related to the CKM matrix. 
The decay modes B — > Z3* =F 7r ± have been proposed 
for use in measurements of sin(2/3 + 7) [2], where (3 = 
arg (— V cd V* b /V td V t * b ) is well measured [3]. In the Stan- 
dard Model the decays B° -> D*-n + and B° -> D*--k+ 
proceed through the b — > cud and b — > ucd amplitudes A c 
and A u . Fig. 1 shows the tree diagrams contributing to 
these decays. The relative weak phase between A u and 
A c in the usual Wolfcnstcin convention [4] is 7. When 
combined with B°B° mixing, this yields a weak phase 
difference of 2/3 + 7 between the interfering amplitudes. 



B° 




D*~ B° 




Here 5* is the strong phase difference between A u and A c , 
and r* = |A U /A C |. Since A u is doubly CKM-suppressed 

v„ h v» 



with respect to A c , one expects r* 



= 0.02. 



We report a study of the CP-violating asymmetry in 
B — > D*^^ decays using the technique of partial re- 
construction, which allows us to achieve a high efficiency 
for the selection of signal events. We use approximately 
twice the integrated luminosity of our previous analysis 
of this process [5], and employ an improved method to 
eliminate a measurement bias, as described in Sec. Ill F 2. 
Many of the tools and procedures used in this analysis 
were validated in a previous analysis dedicated to the 
measurement of the B° lifetime [6] . 

In this analysis, terms of order r* 2 , to which we cur- 
rently have no sensitivity, have been neglected. The in- 
terpretation of the measured asymmetries in terms of 
sin(2/3 + 7) requires an assumption regarding the value 
of r* discussed in Sec. VI. 



II. THE BABAR DETECTOR AND DATASET 



FIG. 1: Feynman diagrams for the Cabibbo-favored decay 
B° —* D*~tt + (left), corresponding to the decay amplitude 
A c , and the Cabibbo-suppressed decay B° — > D*~tt + (right), 
whose amplitude 

In T(4S) — > BB decays, the decay rate distribution for 
B -> D**** is 

e -|A*l/T 

V^(At) = — x [l =p S c sin(AmAf) 

TvC coa(AmAt)} , (1) 

where r is the B° lifetime averaged over the two mass 
eigenstates, Am is the B°B° mixing frequency, and At 
is the difference between the time of the B — > D* =F 7r ± 
(B lec ) decay and the decay of the other B (Stag) in the 
event. The upper (lower) signs in Eq. (1) indicate the 
flavor of the B teig as a B° (B°), while r\ = +1 (—1) and 
C = + (-) for the B ICC final state D*~tt + (D* + tt~). The 
parameters C and S are given by 

C^TT^J' S ± ^-^sin(2/3 + 7± ( T). (2) 
1 + r* z 1 + r* z 



*Also with Universita della Basilicata, Potenza, Italy 
1 Deceased 



The data used in this analysis were recorded with the 
BABAR detector at the PEP-II asymmetric-energy stor- 
age rings, and consist of 211 fb _1 collected on the T(4S) 
resonance (on-resonance sample), and 21 fb" 1 collected 
at an e + e~ center-of-mass (CM) energy approximately 
40 MeV below the resonance peak (off-resonance sam- 
ple). Samples of Monte Carlo (MC) [7] events with 
an equivalent luminosity approximately four times larger 
than the data sample were analyzed using the same re- 
construction and analysis procedure. 

The BABAR detector is described in detail in Ref. [8]. 
We provide a brief description of the main compo- 
nents and their use in this analysis. Charged-particle 
trajectories are measured by a combination of a five- 
layer silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH) in a 1.5-T solenoidal magnetic field. 
Tracks with low transverse momentum can be recon- 
structed in the SVT alone, thus extending the charged- 
particle detection down to transverse momenta of about 
50 MeV/c. We use a ring-imaging Cherenkov detector 
(DIRC) for charged-particle identification and augment it 
with energy-loss measurements from the SVT and DCH. 
Photons and electrons are detected in a CsI(Tl) electro- 
magnetic calorimeter (EMC), with photon-energy reso- 
lution a E /E = 0.023(S/GeV)" 1/4 0.014. The instru- 
mented flux return (IFR) is equipped with resistive plate 
chambers to identify muons. 
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III. ANALYSIS METHOD 

A. Partial Reconstruction of B — > D*^n ± 

In the partial reconstruction of a B — > Z?* =F 7r ± candi- 
date (.B rec ), only the hard (high-momentum) pion track 
TTh from the i? decay and the soft (low-momentum) pion 
track 7r s from the decay D*~ — > £> 7r,T are used. The 
cosine of the angle between the momenta of the B and 
the hard pion in the CM frame is then computed: 



cos 9 Bh 



M 2 B o 



Ml + E CM E h 



2PB\Ph\ 



(3) 



where M x is the nominal mass of particle x [9] , Eh and 
Ph are the measured CM energy and momentum of the 
hard pion, Eqm is the total CM energy of the incom- 

bcams, and p B = 



ing e + e 



M% . 



Events 



are required to be in the physical region \co$,9 Bh \ < 1. 
Given cos 9 Bh and the measured momenta of the Hh and 
ir s , the B four- momentum can be calculated up to an 
unknown azimuthal angle 4> around pu- For every value 
of <f), the expected D four- momentum pd(<I>) is deter- 
mined from four-momentum conservation, and the corre- 
sponding ^-dependent invariant mass m((f>) = \J\pd{4>)\ 2 
is calculated. We define the missing mass m m iss = 
\ [n^max + Wmin] 7 where m max and m m ; n are the max- 
imum and minimum values of m((j>). In signal events, 
TO m j ss peaks at the nominal D° mass M B o, with a gaus- 
sian width of about 3 McV/c 2 (Fig. 2). The m m i ss dis- 
tribution for combinatoric background events is signif- 
icantly broader, making the missing mass the primary 
variable for distinguishing signal from background. The 
discrimination between signal and background provided 
by the m m ; ss distribution is independent of the choice 
of the value of </>. With the arbitrary choice <fr = 0, we 
use four-momentum conservation to calculate the CM D 
and B momentum vectors, which are used as described 
below. 



B. Backgrounds 

In addition to B — > D* T n ± events, the selected event 
sample contains the following kinds of events: 

• B -> D*^p ± . 

• Peaking BB background, defined as decays other 
than B — > D* T p ± , in which the 7r^ and 7r s originate 
from the same B meson, with the 7T S originating 
from a charged D* decay. The m m i ss distribution 
of these events peaks broadly under the signal peak. 

• Combinatoric BB background, defined as all re- 
maining BB background events. 



• Continuum e+e 
s, or c quark. 



qq, where q represents a u, d, 



C. Event Selection 

To suppress the continuum background, we select 
events in which the ratio of the 2nd to the 0th Fox- 
Wolfram moment [10], computed using all charged parti- 
cles and EMC clusters not matched to tracks, is smaller 
than 0.40. Hard-pion candidates are required to be re- 
constructed with at least twelve DCH hits. Kaons and 
leptons are rejected from the tth candidate lists based 
on information from the IFR and DIRC, energy loss in 
the SVT and DCH, or the ratio of the candidate's EMC 
energy deposition to its momentum (E/p). 

We define the D* helicity angle 9 B * to be the angle 
between the flight directions of the D and the B in the 
D* rest frame. Taking advantage of the longitudinal po- 
larization in signal events, we suppress background by 
requiring | cos#£>* | to be larger than 0.4. 

All candidates are required to satisfy rn m ; ss > 
1.81 GcV/c 2 . Multiple candidates are found in 5% of 
the events. In these instances, only the candidate with 
the rrimiss value closest to M D o is used. 



D. Fisher Discriminant 

To further discriminate against continuum events, we 
combine fifteen event-shape variables into a Fisher dis- 
criminant [11] F. Discrimination originates from the fact 
that qq events tend to be jet-like, whereas BB events 
have a more spherical energy distribution. Rather than 
applying requirements to the variable F, we maximize the 
sensitivity by using it in the fits described below. The fif- 
teen variables are calculated using two sets of particles. 
Set 1 includes all tracks and EMC clusters, excluding 
the hard and soft pion candidates; Set 2 is composed of 
Set 1, excluding all tracks and clusters with CM momen- 
tum within 1.25 radian of the CM momentum of the D. 
The variables, all calculated in the CM frame, are 1) the 
scalar sum of the momenta of all Set 1 tracks and EMC 
clusters in nine 20° angular bins centered about the hard 
pion direction; 2) the value of the sphericity, computed 
with Set 1; 3) the angle with respect to the hard pion of 
the sphericity axis, computed with Set 2; 4) the direction 
of the particle of highest energy in Set 2 with respect to 
the hard pion; 5) the absolute value of the vector sum of 
the momenta of all the particles in Set 2; 6) the momen- 
tum \ph\ of the hard pion and its polar angle. 



E. Decay Time Measurement and Flavor Tagging 

To perform this analysis, At and the flavor of the B ta , B 
must be determined. We tag the flavor of the B tag using 
lepton or kaon candidates. The lepton CM momentum 
is required to be greater than 1.1 GeV/c to suppress lep- 
tons that originate from charm decays. If several flavor- 
tagging tracks are present in either the lepton or kaon 
tagging category, the only track of that category used for 
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tagging is the one with the largest value of 9t, the CM 
angle between the track momentum and the momentum 
of the "missing" (unreconstructed) D. The tagging track 
must satisfy cos8t < Ct, where Ct = 0.75 (Ct = 0.50) 
for leptons (kaons), to minimize the impact of tracks orig- 
inating from the decay of the missing D. If both a lepton 
and a kaon satisfy this requirement, the event is tagged 
with the lepton. 

We measure At using At — (z rcc — zt ag )/(7/3c), where 
z ICC (ztag) is the decay position of the B ICC (Stag) along 
the beam axis (z) in the laboratory frame, and the e + e~ 
boost parameter 7/3 is calculated from the measured 
beam energies. To find z lec , we use the ir^ track param- 
eters and errors, and the measured beam-spot position 
and size in the plane perpendicular to the beams (the 
x — y plane). We find the position of the point in space 
for which the sum of the x 2 contributions from the Hh 
track and the beam spot is a minimum. The z coordi- 
nate of this point determines z lcc . The beam spot has 
an r.m.s. size of approximately 120 /im in the horizontal 
dimension (x), 5 /im in the vertical dimension (y), and 
8.5 mm along the beams (z). The average B flight in the 
x — y plane is 30 /im. To account for the B flight in the 
beam-spot-constrained vertex fit, 30 /im are added to the 
effective x and y sizes for the purpose of conducting this 
fit. 

In lepton-tagged events, the same procedure, with the 
7T/j track replaced by the tagging lepton, is used to deter- 
mine z tag - 

In kaon-tagged events, we obtain z tag from a beam- 
spot-constrained vertex fit of all tracks in the event, ex- 
cluding Wh, tt s and all tracks within 1 radian of the D 
momentum in the CM frame. If the contribution of any 
track to the \ 2 °f the vertex is more than 6, the track 
is removed and the fit is repeated until no track fails the 
X 2 < 6 requirement. 

The At error a At is calculated from the results of the 
z rcc and z tag vertex fits. We require |At| < 15 ps and 
a At < 2 ps. 



is the PDF for these events. The PDF for off-resonance 
data is V qq . The parameter values for Vi are different for 
each event type, unless indicated otherwise. Each Vi is 
a product, 

Vi = Mi(m miss ) ^i{F)T/(At,a At ,s t ,s m ), (5) 
where the factors in Eq. (5) are described below. 



1. 



and F PDFs 



The m miss PDF for each event type i is the sum of a 
bifurcated Gaussian plus an ARGUS function [12]: 

Mi(m miss ) = ff Qi{m miss ) + (1 - ff )A i {m misa ), (6) 

where ff is the fractional area of the bifurcated Gaussian 
function. The functions Qi and Ai are 



'exp [-(m - M l ) 2 /2a 2 Ll ] , m < M t 
exp [-(m - M t ) 2 /2a m ] , m > M, ' 



Qi{m) oc 



A(m) oc my 1 - (m/Mf) x 



exp 



(l-{ m /M*) 2 )]6(Mf-m), 



(7) 



(8) 



where Mj is the peak of the bifurcated Gaussian, an 
and am are its left and right widths, ej is the ARGUS 
exponent, Mf is its end point, and 6 is the step func- 
tion. The proportionality constants are such that each 
of these functions is normalized to unit area within the 
TO m j ss range. The m m i ss PDF of each event type has dif- 
ferent parameter values. 

The Fisher discriminant PDF Ti for each event type is 
parameterized as the sum of two Gaussians. The parame- 
ter values of Tw-rx-, ^D'p, -^poak, and JFcomb are identical. 



2. Signal At PDFs 



F. Probability Density Function 

The probability density function (PDF) depends on 
the variables m m i ss , At, a At, F, s t , and s m , where s t = 
1 (—1) when the B tag is identified as a B° (B°), and 
s m = 1 (—1) for "unmixed" ("mixed") events. An event 
is labeled unmixed if the 717, is a ir~(n + ) and the B tag is 
a B°(B°), and mixed otherwise. 

The PDF for on-resonance data is a sum over the PDFs 
of the different event types: 



(4) 



where the index i — {D*ir, D*p, peak, comb, qq} indicates 
one of the event types described above, fi is the relative 
fraction of events of type i in the data sample, and Vi 



The At PDF T^. n (At, a A t, s t , s m ) for signal events 
corresponds to Eq. 1 with 0{r* 2 ) terms neglected, modi- 
fied to account for several experimental effects, described 
below. 

The first effect has to do with the origin of the tagging 
track. In some of the events, the tagging track originates 
from the decay of the missing D. These events are labeled 
"missing-D tags" and do not provide any information re- 
garding the flavor of the -Btag- In lepton-tagged events, 
we further distinguish between "direct" tags, in which 
the tagging lepton originates directly from the decay of 
the Stag, and "cascade" tags, where the tagging lepton 
is a daughter of a charmed particle produced in the -Btag 
decay. Due to the different physical origin of the tagging 
track in cascade and direct tags, these two event cate- 
gories have different mistag probabilities, defined as the 
probability to deduce the wrong B flavor from the charge 
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of the tagging track. In addition, the measured value of 
z tag m cascade-lcpton tags is systematically larger than 
the true value, due to the finite lifetime of the charmed 
particle and the boosted CM frame. This creates a corre- 
lation between the tag and vertex measurements that wc 
address by considering cascade-lepton tags separately in 
the PDF. In our previous analysis [5] we corrected for the 
bias of the S' ± parameters caused by this effect and in- 
cluded a systematic error due to its uncertainty. In kaon 
tags, z tag is determined using all available B tag tracks, so 
the effect of the tagging track on the z tag measurement 
is small. Therefore, the overall bias induced by cascade- 
kaon tags is small, and there is no need to distinguish 
them in the PDF. 

The second experimental effect is the finite detector 
resolution in the measurement of At. We address this by 
convoluting the distribution of the true decay time dif- 
ference At tr with a detector resolution function. Putting 
these two effects together, the At PDF of signal events 
is 

T^(Ai,<7At,s t ,0 = + AfD'^EA', x 

3 

J dAt tI Tg.^Aitr, s t , sj Tl j D , n {At - Attr, <7 At ),(9) 

where Ae^.^ is half the relative difference between the 
detection efficiencies of positive and negative leptons or 
kaons, the index j = {dir, cas, miss} indicates di- 
rect, cascade, and missing- D tags, and f° D » n is the frac- 
tion of signal events of tag-type j in the sample. Wc 
set = 1 - f c D % - for lepton tags, with 

the value = 0.12 ± 0.02 obtained from the MC 

simulation. For kaon tags /f, 1 ^ = 0. The function 
T^.^Aitr, s t , s m ) is the At tr distribution of tag-type j 
events, and IZ^.^At— Ai tr , a At) is their resolution func- 
tion, which parameterizes both the finite detector reso- 
lution and systematic offsets in the measurement of Az, 
such as those due to the origin of the tagging particle. 
The parameterization of the resolution function is de- 
scribed in Sec. IIIF4. 

The functional form of the direct and cascade tag Ai tr 
PDFs is 

e -\At t r\/T D ,„ 

T^JAt tr ,s t ,sJ = — x 

|l-s t AJ D , w 

+s m (1 - 2uj J DtJT ) cos(A™_ D « 7r At t r) 
-S j D , v sin(A™ I) « Tr At tr )|, (10) 

where j = {dir, cas}, the mistag rate u^.^ is the proba- 
bility to misidentify the flavor of the B tag averaged over 
B° and _B°, and Acj : ' D , 7r is the B° mistag rate minus the 
B° mistag rate. The factor S J D ,^ describes the effect of 
interference between b — > ucd and b — > cud amplitudes in 



both the B rec and the B tag decays: 

^d*tt = (l-2^, 7r )(s t a D , 7r + s m c D , 7r ) 

+s i sJ> D . 7I {l-s i Auj j D ,„), (11) 

where cid*it, ^d*tt, and Cd*k are related to the physical 
parameters through 

aD'n = 2r* sin(2/3 + 7) cos<5*, 

b D , n = 2r'sin(2/3 + 7)cos(5', 

c D , n = 2cos(2/3 + 7)(r*sin(5* -r'sinS'), (12) 

and r' (5') is the effective magnitude of the ratio (effec- 
tive strong phase difference) between the b — > ucd and 
b — > cud amplitudes in the B tag decay. This parameteri- 
zation is good to first order in r* and r' . In the following 
we will refer to the parameters cid'tt, bw-n, cd'k and 
related parameters for the background PDF as the weak 
phase parameters. Only aw-n and bo*-* are related to CP 
violation, while cx>. T can be non-zero even in the absence 
of CP violation when 2/3 + 7 = 0. The inclusion of r' and 
5' in the formalism accounts for cases where the B tag un- 
dergoes a b — ► ucd decay, and the kaon produced in the 
subsequent charm decay is used for tagging [13]. We ex- 
pect r' ~ 0.02. In lepton-tagged events r' = (and hence 
bD»TT = 0) because most of the tagging leptons come from 
B semileptonic decays to which no suppressed amplitude 
with a different weak phase can contribute. 
The At tr PDF for missing- D tags is 

T™(At tI ,s t ,sJ = e "'^°" (l + am (l-2p D . ir ) 

- 2s t s n b D *TT sin(AmD«7rAt t r)|,(13) 

where pd'-k is the probability that the charge of the tag- 
ging track is such that it results in a mixed flavor mea- 
surement. In this analysis, we have neglected the term 
proportional to sin(Am_D« 7r Ai tr ) of Eq. 13. The system- 
atic error on bo*-n due to this approximation is negligible 
due to the small value of /j£J^ reported below. 



3. Background At PDFs 

The At PDF of B -> D*^p ± has the same functional 
form and parameter values as the signal PDF, except that 
the weak phase parameters cld* p , br)* P , and cd* P are set 
to and are later varied to evaluate systematic uncer- 
tainties. The validity of the use of the same parameters 
for 7^)* p and Tp* n is established using simulated events, 
and stems from the fact that the TTh momentum spectrum 
in the B — ► Z?* =F p ± events that pass our selection criteria 
is almost identical to the signal spectrum. 

The At PDF of the peaking background accounts sep- 
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arately for charged and neutral B decays: 

T p ' cak (Ai, (J At: s t ,sj = (1 + s t Ae pcak ) {r° poak 
+ J dAt tT T+ eajk (At tT ,s t ,sJ x 

ft+ak( A *- A *tr>*At)}> 



(14) 



where T° pcak has the functional form of Eq. (9) and 
the subsequent expressions, Eqs. (13-12), but with all 
Z?*7r-subscripted parameters replaced with their peak- 
subscripted counterparts. The integral in Eq. (14) ac- 
counts for the contribution of charged B decays to the 
peaking background, with 



e -|Ai tr |/T+ ak 



peak / \ 

+ (l-« t A< eak ), (15) 



pcakv »*>-w , 

peak 

and 7?. pcak (At — At tI , <JAt) being the three-Gaussian res- 
olution function for these events described below. 

The Combinatoric BB background PDF 7^,' omb is sim- 
ilar to the signal PDF, with one substantial difference. 
Instead of parameterizing T^ olnh with the four parame- 
ters /mmb; ^comb- Aw comb> Pcomb, we use the set of three 
parameters 

i^dir 

I _ dir /-, _ rdir \ . ./ comb 

^comb ^comb V ± J comb ) ' r> > 



A^comb — Aw comb (1 - /comb)) 
ftomb = ,/comb(l — 2 Pcomb)- 



(16) 



With these parameters and /co mb = the combinatoric 
BB background At PDF becomes 

T c'omb( A *> "At, s t ,sj = (1 + s t Ae comb ) X 

dAt tI T comh (At tI , s t , s m ) TZ comh (At - At tr , a At) ,(17) 



where 7£ C omb(At — At tr , &At) is the 3-Gaussian resolution 
function and 



% omb (At t r, S t , S m ) = 



-|At tr |/r comb 
4T CO mb 



1 - s t Auil 



comb 



+s m ^comb + s m (1 - 2uj' comh ) cos(Am comb Ai t r) 
-5 comb sin(Am comb Ai t r) \, (18) 



with 

<5 C omb = (1 — 2tJ comb ) (s t a comb + S m C com b) 

+ s t s m 6 comb (l - s t Au}' comh ). (19) 

As in the case of Tb* p , the weak phase parameters of 
the peaking and combinatoric background (a pe ak, &peak, 
Cpcak and a comb , 6 com b, c comb ) are set to and are later 
varied to evaluate systematic uncertainties. Parameters 
labeled with superscripts "peak" or "comb" are empirical 



and thus do not necessarily correspond to physical pa- 
rameters. In general, their values may be different from 
those of the £>*7r-labeled parameters. 

The PDF T q q for the continuum background is the sum 
of two components, one with a finite lifetime and one with 
zero lifetime: 

Tqq{^t,o&t,s t ) = (l + s t Ae q q) J dAt tT T q g(At tl , s t , S m ) 
x 1l q q{At- At tT ,a At ), (20) 

with 

e -|At tr |/r,- 



Tqq(At tl ,S t ) = (1-.4) 



(1 - s t Aw, 



qq, , V- "t qq) 

VTqq 

+ fqq6(At tl ), (21) 



where ff l7! is the fraction of zero-lifetime events. 
•' qq 



4- Resolution Function Parameterization 

The resolution function for events of type i and op- 
tional secondary type j (j = {dir, cas, miss} for lepton- 
tagged signal events and j = {+, 0} for the peaking and 
combinatoric BB background types) is parameterized as 
the sum of three Gaussians: 

ni(t r ,<TAt) = f n iQ n i(t r ,CT At) 

+ nG i(t r ,a At ), (22) 

where t r = At — At tr is the residual of the At measure- 
ment, and Q n j, G w \-, and Q°\ are the "narrow", "wide", 
and "outlier" Gaussians. The narrow and wide Gaussians 
have the form 



G k i(t r ,a A t) 



2TTS k j a At 



cxp - 



(tr 



b k3 



2( S fc ^A t ) 2 



, (23) 



where the index k takes the values k = n,w for the nar- 
row and wide Gaussians, and b k \ and s k \ are parameters 
determined by fits, as described in Sec. Ill G. The outlier 
Gaussian has the form 



G°\{t r ,OAt) 



1 



2irs oJ - 



exp 



2(s1) 2 



(24) 



where in all nominal fits the values of b°\ and s°\ are 
fixed to ps and 8 ps, respectively, and are later varied 
to evaluate systematic errors. 
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G. Analysis Procedure 



IV. RESULTS 



The analysis is carried out with a series of unbinned 
maximum-likelihood fits, performed simultaneously on 
the on- and off-resonance data samples and indepen- 
dently for the lepton-tagged and kaon-tagged events. The 
analysis proceeds in four steps: 

1. In the first step, we determine the parameters 
Id*p + /d*tt, /peak, and / comb of Eq. (4). In or- 
der to reduce the reliance on the simulation, we 

also obtain in the same fit the parameters of 
Eq. (6), e q q of Eq. (8), ol for the signal m m ; ss PDF 
(Eq. (7)), and all the parameters of the Fisher dis- 
criminant PDFs. This is done by fitting the data 
with the PDF 



Vi = Mi(m miss )Fi(F), 



(25) 



instead of Eq. (5); i.e. by ignoring the time de- 
pendence. The fraction f qq of continuum events is 
determined from the off-resonance sample and its 
integrated luminosity relative to the on-resonance 
sample. All other parameters of the A4i PDFs and 
the value of /r>*7r/(/i)*7r + Jd* p ) = 0.87 ± 0.03 are 
obtained from the MC simulation. 

2. In the second step, we repeat the fit of the first 
step for data events with cos 6t > Ct, to ob- 
tain the fraction of signal events in that sample. 
Given this fraction and the relative efficiencies for 
direct, cascade, and missing- D signal events to sat- 
isfy the cos 9t < Ct requirement, we calculate 
/jSi™ = 0.011 ± 0.001 for lepton-tagged events and 
f%?* = 0.055 ± 0.001 for kaon-tagged events. We 
also calculate the value of pd** from the frac- 
tions of mixed and unmixed signal events in the 
cos 9t > Ct sample relative to the cos 9t < Ct 
sample. 

3. In the third step, we fit the data events in the 
sideband 1.81 < m miss < 1.84 GeV/c 2 with the 
3-dimensional PDFs of Eq. (5). The parameters 
of A4i(m m i ss ) and Fi{F), and the fractions /, are 
fixed to the values obtained in the first step. From 
this fit we obtain the parameters of ^.' omb , as well 
as those of T^~. 

4. In the fourth step, we fix all the parameter values 
obtained in the previous steps and fit the events in 
the signal region m m ; ss > 1.845 GeV/c 2 , determin- 
ing the parameters of 77-,, ^ and 77^. Simulation 
studies show that the parameters of 7^' omb are in- 
dependent of m m i ss , enabling us to obtain them in 
the sideband fit (step 3) and then use them in the 
signal-region fit. The same is not true of the 77= 
parameters; hence they are free parameters in the 
signal-region fit of the last step. The parameters of 
7^'cak are obtained from the MC simulation. 



The fit of step 1 finds 18710 ± 270 signal B -> D* t tt ± 
events in the lepton-tag category and 70580 ± 660 in the 
kaon-tag category. The m m ; ss and F distributions for 
data are shown in Figs. 2 and 3, with the PDFs overlaid. 




miss 



1.87 l.S 

(GeV/c 2 ) 



FIG. 2: The m m i ss distributions for on-resonance lepton- 
tagged (top) and kaon-tagged (bottom) data. The curves 
show, from bottom to top, the cumulative contributions of 
the continuum, peaking BB, combinatoric BB, B — » D* 
and B — > D*^ 1 ^ PDF components. 



The results of the signal region fit (fourth step) are 
summarized in Table I, and the plots of the At distri- 
butions for the data are shown in Fig. 4 for the lepton- 
tagged and the kaon-tagged events. The goodness of the 
fit has been verified with the Kolmogorov-Smirnov test 
and by comparing the likelihood obtained in the fit with 
the likelihood distribution of many parameterized MC 
experiments generated with the PDF's obtained in the 
fit on the data. Fig. 5 shows the raw, time-dependent 
CP asymmetry 



A(At) 



N St=1 (At) - N^^jAt) 
N St=1 {At)+N s -^(At)- 



(26) 



In the absence of background and with high statistics, 
perfect tagging, and perfect At measurement, A(At) 
would be a sinusoidal oscillation with amplitude o^.j. 
For presentation purposes, the requirements m m j ss > 
1.855 GeV/c 2 and F < were applied to the data plot- 
ted in Figs. 4 and 5, in order to reduce the background. 
These requirements were not applied to the fit sample, 
so they do not affect our results. 

The fitted values of Am reported in Table I are in good 
agreement with the world average (0.502 ± 0.007) ps _1 




F 

FIG. 3: The F distributions for on-resonance lepton-tagged 
(top) and kaon-tagged (bottom) data. The contributions of 
the BB (dashed-dotted line) and the continuum (dashed line) 
PDF components are overlaid, peaking at approximately —0.6 
and —0.1, respectively. The total PDF is also overlaid. 

[9]. The fitted values of the B° lifetime need to be 
corrected for a bias observed in the simulated sam- 
ples, At = Tfn — r gen = (—0.03 ± 0.02) ps for the 
lepton-tag and At = (—0.04 ± 0.02) ps for the kaon- 
tag events. After this correction, the measured life- 
times, t(B°) = (1.48 ± 0.02 ± 0.02) ps and t(B°) = 
(1.49 ± 0.01 ± 0.04) ps for the lepton-tag and kaon-tag, 
respectively, are in reasonable agreement with the world 
average t(B°) = (1.536 ± 0.014) ps [9]. The correla- 
tion coefficients of a D , v (<^>»„) with Am and t(B°) are 
-0.021 and 0.019 (-0.060 and -0.056). 



V. SYSTEMATIC STUDIES 

The systematic errors are summarized in Table II. 
Each item below corresponds to the item with the same 
number in Table II. 

1. The statistical errors from the fit in Step 1 are prop- 
agated to the final fit. This also includes the sys- 
tematic errors due to possible differences between 
the PDF line shape and the data points. 

2. The statistical errors from the m m i ss sideband fit 
(Step 3) are propagated to the final fit (Step 4). 

3-4. The statistical errors from the Step 2 fits are prop- 
agated to the final fit. 
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A t (ps) 



FIG. 4: At distributions for the lepton-tagged (a-d) and kaon- 
tagged (e-h) events separated according to the tagged flavor 
of Stag and whether they were found to be mixed or unmixed: 
a,e) B° unmixed, b,f) B° unmixed, c,g) B° mixed, d,h) B° 
mixed. The solid curves show the PDF, calculated with the 
parameters obtained by the fit. The PDF for the total back- 
ground is shown by the dashed curves. 



5. The statistical errors associated with the parame- 
ters obtained from MC are propagated to the fi- 
nal fit. In addition, the full analysis has been per- 
formed on a simulated sample to check for a pos- 
sible bias in the weak phase parameters measured. 
No statistically significant bias has been found and 
the statistical uncertainty of this test has been as- 
signed as a systematical error. 

6. The effect of uncertainties in the beam-spot size on 
the vertex constraint is estimated by increasing the 
beam spot size by 50 fan. 

7. The effect of the uncertainty in the measured length 
of the detector in the z direction is evaluated by 
applying a 0.6% variation to the measured values 
of At and a At- 
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TABLE I: Results of the fit to the lepton- and kaon-tagged events in the signal region 1.845 < m m i ss < 1.880 GeV/c 2 . Errors 
are statistical only. See Sections IIIF2, IIIF3, and IIIF4 for the definitions of the symbols used in this table. 





Lepton tags 


Kaon tags 


Parameter description 


Parameter 


Value 


Parameter 


Value 


Signal weak phase par. 


On,. 


-0.042 ± 0.019 


K 

7T 


-0.025 ± 0.020 
-0.004 ±0.010 




CjJ * Tj- 


-0.019 ±0.022 




-0.003 ±0.020 


Signal At PDF 




0.518 ± 0.010 ps _i 




0.4911 ± 0.0076 ps -1 




T D*7T 


1.450 ±0.017 ps 


T~D*n 


1.449 ±0.011 ps 




, ,dir 


0.010 ± 0.006 




0.2302 ± 0.0035 
-0.0181 ±0.0068 






0.027 ±0.010 


Ae]j*Tr 


-0.0070 ± 0.0073 


Signal resolution function 


r n cas 
°D*ir 
rw cas 
°D*ir 


-0.58 ±0.16 
0.23 ±2.01 








rn dir 

°D*7T 


0. (fixed) 


D'TT 


-0.255 ±0.013 




i w dir 

U D*7T 


0. (fixed) 


h w 


-2.07 ±0.48 




rn dir 

J D * 7T 


0.978 ± 0.008 


rn 

J D * 7T 


0.969 ± 0.007 




dir 


0. (fixed) 


fr>* 


0.000 ± 0.001 




n dir 


1.080 ± 0.033 


s" 


1.029 ± 0.023 




iti dir 

Z3 * 7T 


5.76 ± 1.44 


b D*TT 


4 35 + 40 


Continuum At PDF 




1.26 ± 0.32 ps 




0.707 ± 0.048 ps 






0.340 ± 0.009 


f 


0.045 ± 0.022 
0.311 ±0.006 




f*_ 
Jqq 


0.815 ± 0.064 


Jqq 


0.820 ±0.015 


Continuum resolution function 


b n - 

u qq 


0.026 ± 0.048 


b n - 

u qq 


0.017 ± 0.005 




b w - 

u qq 


-0.39± 0.23 


bqq 


-0.043 ± 0.043 




rn 
J qq 


0.65 ± 0.12 


rn 
Jqq 


0.858 ± 0.014 




f°- 
J qq 


0.068 ± 0.014 


f°- 
J qq 


0.018 ± 0.001 




s n - 

*qq 


0.929 ± 0.078 


"qq 


1.064 ± 0.008 




b iq 


1.81 ±0.28 


»qq 


2.267 ± 0.099 



8. To evaluate the effect of possible misalignments in 
the SVT, signal MC events are reconstructed with 
different alignment parameters, and the analysis is 
repeated. 

9-11. The weak phase parameters of the B — > _D* T p ± , 
peaking, and combinatoric BB background are 
fixed to in the fits. To study the effect of pos- 
sible interference between b — > ucd and b — > cud 
amplitudes in these backgrounds, their weak phase 
parameters are varied in the range ±0.04 and the 
Step-4 fit is repeated. We take the largest variation 
in each weak phase parameter as its systematic er- 
ror. 

12. In the final fit, we take the values of the parame- 
ters of T^ eak from a fit to simulated peaking BB 
background events. The uncertainty due to this is 
evaluated by fitting the simulated sample, setting 
the parameters of T' k to be identical to those of 



13. The uncertainty due to possible differences between 
the At distributions for the combinatoric back- 
ground in the m m i ss sideband and signal region is 
evaluated by comparing the results of fitting the 



simulated sample with the ^.' omb parameters taken 
from the sideband or the signal region. 

14. The ratio /d*p//d*tt is varied by the uncertainty 
in the corresponding ratio of branching fractions, 
obtained from Ref. [9]. 



VI. PHYSICS RESULTS 

Summarizing the values and uncertainties of the weak 
phase parameters, we obtain the following results from 
the lepton-tagged sample: 

a l D ^ = -0.042 ±0.019 ±0.010, 

c e D , v = -0.019 ±0.022 ±0.013. (27) 

The results from the kaon-tagged sample fits are 

a%, n = -0.025 ± 0.020 ±0.013, 
b%, n = -0.004 ±0.010 ±0.010, 
c%, n = -0.003 ± 0.020 ±0.015. (28) 
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TABLE II: Systematic errors in </ D *^ and c l D *^ for lepton-tagged events and a^,, n , bp* n , and Cq*^ for kaon-tagged events. 



Source 


Error (xl0~ 2 ) 


Lepton tags 


Kaon tags 


i i 


n K h K r K 

LJ TV LJ TV U TV 


1. Step 1 fit 

2. Sideband statistics 

o rmiss 

JD»7T 

4. pc* w 


0.04 0.04 
0.08 0.08 
0.02 0.02 
0.02 0.02 


0.10 0.04 0.04 
0.40 0.12 0.44 
0.02 negl. negl. 
0.02 negl. negl. 


5. MC statistics 


0.60 0.82 


0.68 0.34 0.70 


6. Beam spot size 

7. Detector z scale 

8. Detector alignment 


0.10 0.10 
0.03 0.03 
0.25 0.55 


0.07 0.13 0.06 
0.02 negl. 0.03 
0.25 0.13 0.41 


9. Combinatoric background weak phase par. 

10. Peaking background weak phase par. 

11. D* p weak phase par. 


0.25 0.22 
0.36 0.38 
0.53 0.52 


0.80 0.56 0.72 
0.29 0.17 0.27 
0.57 0.58 0.58 


12. Peaking background 

13. Signal region/sideband difference 

14. B(B D^p*) 


0.21 0.31 
negl. negl. 
0.17 0.33 


0.21 0.41 0.31 
0.04 0.03 0.05 
0.17 0.22 0.33 


Total systematic error 


1.0 1.3 


1.4 1.0 1.5 


Statistical uncertainty 


1.9 2.2 


2.0 1.0 2.0 



Combining the results for lepton and kaon tags gives the 
amplitude of the time-dependent CP asymmetry, 

ao*TT = 2r* sin(2/3 + 7) cos 5* 

= -0.034 ±0.014 ± 0.009, (29) 

where the first error is statistical and the second is sys- 
tematic. The systematic error takes into account cor- 
relations between the results of the lepton- and kaon- 
tagged samples coming from the systematic uncertain- 
ties related to detector effects, to interference between 
b — > ucd and b — > cud amplitudes in the backgrounds 
and from B(B — > £>* =F p ± ). This value of an*-* deviates 
from zero by 2.0 standard deviations. 

Previous results of time-dependent CP asymmetries re- 
lated to 2/3 ±7 appear in Ref. [5, 14]. This measurement 
supersedes the results of the partial reconstruction anal- 
ysis reported in Ref. [5] and improves the precision on 
an* n an d Cd*-k with respect to the average of the pub- 
lished results. 

We use a frequentist method, inspired by Ref. [15], to 
set a constraint on 2/3 + 7. To do this, we need a value 
for the ratio r* of the two interfering amplitudes. This is 
done with two different approaches. 

In the first approach, to avoid any assumptions on the 
value of r* , we obtain the lower limit on | sin(2/3 + 7)! as 
a function of r* . 

We define a x 2 function that depends on r* , 2/3 + 7, 
and 6* : 

3 

X 2 (r*,2f3 + 1 ,5*) = ]T Ax^' 1 Ax k , (30) 

j;k = l 

where Axj is the difference between the result of our mea- 
surement of a^. n , a e D , n , or c l D ,^ (Eqs. (28) and (27)) 
and the corresponding theoretical expressions given by 



Eq. (12). We fix r* to a trial value r°. The measure- 
ments of b 1 ^,^ and c^,^ arc not used in the fit, since 
they depend on the unknown values of r' and 5' . The 
measurement error matrix V is nearly diagonal, and ac- 
counts for correlations between the measurements due to 
correlated statistical and systematic uncertainties. We 
minimize % 2 as a function of 2/3 + 7 and 5* , and obtain 
Xmiw t ne minimum value of x 2 . 

In order to compute the confidence level for a given 
value x of 2/3 + 7, we perform the following procedure: 

1. We fix the value of 2/3 + 7 to x and minimize % 2 
as a function of 5*. We define Xmin( x ) to be the 
minimum value of the x 2 in this fit, and Sf Qy to 
be the fitted value of 8*. We define Ax 2 (x) = 

Xmin Xmin ' 

2. We generate many parameterized MC experiments 
with the same sensitivity as the data sample, taking 
into account correlations between the observables, 
expressed in the error matrix V of Eq. (30). To 
generate the observables a^, n , a e D »^, and Cjj,^, we 
use the values (2/3 + 7) = x, r* = r° and S* = 
5t oy . For each experiment we calculate the value of 
Ax 2 (x), computed with the same procedure used 
for the experimental data. 

3. We interpret the fraction of these experiments for 
which Ax 2 (x) is smaller than Ax 2 (x) in the data 
to be the confidence level (CL) of the lower limit 
on (2/3 + 7) = x. 

The resulting 90% CL lower limit on | sin(2/3 ±7)! as 
a function of r* is shown in Fig. 6. The x 2 function is 
invariant under the transformation 2/3 + 7 — * tt/2 + 5* 
and 5* — > ir/2 — 2/3 + 7. The limit shown in Fig. 6 is 
always the weaker of these two possibilities. 
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FIG. 5: Raw asymmetry for (a) lepton-tagged and (b) kaon- 
tagged events. The curves represent the projections of the 
PDF for the raw asymmetry. A nonzero value of ao'-n would 
show up as a sinusoidal asymmetry, up to resolution and back- 
ground effects. The offset from the horizontal axis is due to 
the nonzero values of Aeo'-w and Aud**- 



In the second approach, we estimate r* as originally 
proposed in Ref. [2], and assume SU(3) flavor symme- 
try. With this assumption, r* can be estimated from 
the Cabibbo angle Oc, the ratio of branching fractions 
B(B° -> D*jn-)/B{B° -> D*-n+) = (5.4±|3 ± 0.7) x 
10 -3 [16], and the ratio of decay constants /d*//d* = 
1.10 ±0.02 [17], 



' B(B a -> D* 



y B{B° D*-n+) f D 
yielding the measured value 



- — tan(6>c) 



0.015 



+0.004 
-0.006- 



(31) 

(32) 
4>n + ), for 



This value depends on the value of B(Df 
which we use our recent measurement [18]. 

Equation (31) has been obtained with two approxima- 
tions. In the first approximation, the exchange diagram 




0.01 



0.02 



0.03 



0.04 



0.05 



FIG. 6: Lower limit on | sin(2/3 + 7)| at 90% CL as a function 
of r*, for r* > 0.001. 



amplitude E contributing to the decay B° — > D* + tt~ 
has been neglected and only the tree-diagram ampli- 
tude T has been considered. Unfortunately, no reli- 
able estimate of the exchange term for these decays ex- 
ists. The only decay mediated by an exchange diagram 
for which the rate has been measured is the Cabibbo- 
allowed decay B° -> DjK + . The average of the BABAR 
and Belle branching fraction measurements [16, 19] is 
(3.8 ± 1.0) x 10~ 5 . This yields the approximate ra- 
tio B{B° -> D-K + )/B(B° -> D-TT+) ~ 10~ 2 , which 
confirms that the exchange diagrams are strongly sup- 
pressed with respect to the tree diagrams. Detailed 
analyses [20] of the B — » Dir and B — > D*ir decays 
in terms of the topological amplitudes conclude that 
\E'/T'\ = 0.12 ±0.02 for B° -> D~n+_ and \E/f_\ < 0.10 
for B° -> D*~tt+ decays, where E' , E and T' , T are the 
exchange and tree amplitudes for these Cabibbo-allowed 
decays. We assume that a similar suppression holds for 
the Cabibbo-suppressed decays considered here. 

The second approximation involves the use of the ratio 
of decay constants /d*//d* to take into account SU(3) 
breaking effects and assumes factorization. We attribute 
a 30% relative error to the theoretical assumptions in- 
volved in obtaining the value of r* of Eq. (32), and use 
it as described below. 

We add to the x 2 of Eq. (30) the term A 2 (r*) that 
takes into account both the Gaussian experimental errors 
of Eq. (32) and the 30% theoretical uncertainty according 
to the prescription of Ref. [21]: 



A 2 (r*) 



- 1.3 r* 



0.004 



r* -0.7r*' 
0.006 



> 0.3 



IM<0.3 , (33) 



< -0.3 , 



where £ r . = (r* - r * meas )/r wea8 . 

To obtain the confidence level we have repeated the 
procedure described above with the following changes. 
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°0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

|sin(2(3+y)l 

FIG. 7: The shaded region denotes the allowed range of 
j sin(2/3 + 7)| for each confidence level. The horizontal lines 
show, from top to bottom, the 68% and 90% CL. 

To compute Xmin we minimize x 2 as a function of 2/3 + 7, 
r* and 5*. The value Xmin( x ) ^ s obtained minimizing x 2 
as a function of r* and S* , having fixed 2/3 ± 7 to a given 
value x. We define 5^ oy and r* toy to be the fitted value of 
6* and r* in this fit. To generate the observables a^,^, 
a>jj* n , and c^,,^ in the parameterized MC experiments, 
we use the values (2/3 + 7) = x, r* = r* toy and 6* = SI . 

The confidence level as a function of | sin(2/3 + 7)) is 
shown in Fig. 7. We set the lower limits | sin(2/3 + 7)! > 
0.62 (0.35) at 68% (90%) CL. The implied probability 
contours for the apex of the unitarity triangle, parame- 
terized in terms of p and fj defined in Ref. [4], appear in 
Fig. 8. 

VII. SUMMARY 

We present a measurement of the time-dependent CP 
asymmetries in a sample of partially reconstructed B° — > 
D* + tt~ events. In particular, we have measured the pa- 
rameters related to 2/3 + 7 to be 

cld'tt = 2r* sin(2/3 + 7) cos 5* 

= -0.034 ±0.014 ±0.009 (34) 

and 

c* D ,„ = 2r*cos(2/3 + 7 )sin<5* 

= -0.019 ±0.022 ±0.013, (35) 

where the first error is statistical and the second is sys- 
tematic. We extract limits as a function of the ratio r* 
of the b — > ucd and b — * cud decay amplitudes. With 
some theoretical assumptions, we interpret our results in 
terms of the lower limits | sin(2/3 + 7)! > 0.62 (0.35) at 
68% (90%) CL. 




FIG. 8: Contours of constant probability (color-coded in per- 
cent) for the position of the apex of the unitary triangle to 
be inside the contour, based on the results of Fig. 7. The 
cross represents the value and errors on the position of the 
apex of the unitarity triangle from the CKMFitter fit using 
the "ICHEP04" results excluding this measurement [22] . 
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